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ABSTRACT 

We study the peculiar velocity field inferred from the Mark III spirals using a new 
method of analysis. We estimate optimal values of TuUy-Fisher scatter and zero- 
point offset, and we derive the 3-dimensional rms peculiar velocity (ctu) of the galax- 
ies in the samples analysed. We check our statistical analysis using mock catalogs 
derived from numerical simulations of CDM models considering measurement un- 
certainties and sampling variations. Our best determination for the observations is 
av = (660 ± 50)km/s. We use the linear theory relation between a^ , the den- 
sity parameter Q, and the galaxy correlation function i^(r) to infer the quantity 
/3 = riP-^ /b = 0.601q :[^ where b is the linear bias parameter of optical galaxies and the 
uncertainties correspond to bootstrap resampling and an estimated cosmic variance 
added in quadrature. Our findings are consistent with the results of cluster abundances 
and redshift space distortion of the two-point correlation function. These statistical 
measurements suggest a low value of the density parameter fi ~ 0.4 if optical galaxies 
are not strongly biased tracers of mass. 

Key words: galaxies — distance scale — velocity field — cosmology: density param- 
eter 



1 INTRODUCTION 



Recent developements on extragalactic distance indicators 
(Djorgovski and Davis, 1987; Dressier et al., 1987a) allow 
us to study the peculiar galaxy velocity field in the local 
Universe up to ~ 50 — 100/i~^Mpc (see Giovanelli, 1997, 
or Strauss and Willick, 1995 for a review). These measure- 
ments of peculiar velocities provide direct probes of the mass 
distribution in the Universe and set constraints on models 
of large-scale structure formation. By comparing the mass 
distribution implied by the velocity field with the observed 
distribution of galaxies the density parameter Q. can be esti- 
mated. Nevertheless, Q, may only be determined within the 
uncertainty of the bias parameter b {b — 1/as is the in- 
verse of the root mean square mass fluctuations in spheres 
of radius= 8Mpc /i"^) through the factor f3 = n°-^/b. 

Bertschinger & Dekel (1989) developed the POTENT 
method whereby the mass distribution may be reconstructed 
by using the analog of the Bernoulli equation for irrota- 
tional flows. This method was used to analyse the peculiar 
velocity fleld out to 60/i~^Mpc (Bertschinger et al., 1989). 
Dekel et al. (1993) compared the previously determined ve- 
locity field with the observed distribution of galaxies con- 
cluding that il^'^/b ~ 1 provides the best-fitting to the 
data. Analysis of the velocity tensor (Gorski, 1988; Groth, 
Juszkiewics and Ostriker, 1989) also provide useful insights 
on the velocity field. Zaroubi et al. (1997) reconstructed the 
large-scale power spectrum from the velocity tensor of the 
Marklll data and found it consistent with a CDM model 
with asQ°-^ ~ 0.8 although a different resuh, agil"-^ ~ 0.35, 
is found by Kashlinsky (1997) in a similar analysis. 



Relations between root mean square mass fluctuation 
in a given scale and density parameter may also be obtained 
in studies of cluster abundances in different cosmological 
models. These analysis place constraints of the form f2'^/6 ~ 
0.4 — 0.6 with a ~ 0.4 — 0.6 in a variety of cold and mixed 
dark matter models (see Eke et al., 1998, and Gross et al., 
1998). Similarly, studies of redshift space distortions of the 
galaxy two point correlation function also provide a useful 
restriction to the parameter /3, as for instance Ratcliffe et 
al. (1997) who find (3 ~ 0.5. 

In this paper we study the peculiar velocity field 
through a statistical analysis of observational data taken 
from the Mark III catalog. Data characteristics are presented 
in section 2. Section 3 provides an outline of the statistical 
procedure. In section 4 we analyse the Mark III spirals, and 
in section 5 we test our procedure using mock catalogs ac- 
cording to different observers in fully non-linear numerical 
simulations. In section 6 we provide a determination of the 
parameter /3. 



2 DATA 

We use samples of spiral galaxies taken from the Mark III 
catalog (Willick et al., 1995; Willick et al., 1996; Willick et 
al., 1997) to analyze the peculiar velocity flow. This Cata- 
log lists TuUy-Fisher and Z)„ — a distances and radial ve- 
locities for spiral, irregular, and elliptical galaxies. For spi- 
ral galaxies, the velocity parameter rj — LogAV — 2.5 is 
determined either from HI profiles or from optical Ha rota- 
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tion curves. The TuUy-Fisher (TF) relations and their corre- 
sponding scatters for the different samples of spiral galaxies 
are given by Willick et al. (1997) and are shown in Table 1, 
where the absolute magnitude M satisfies M = m — Slogcz. 
The galcixy apparent magnitudes m of the TuUy-Fisher dis- 
tances are corrected for Galactic extinction, inclination and 
redshift (see Wilhck et al. 1997 for details). 

The selection bias in the calibration of the forward 
TF relation can be corrected once the selection function is 
known. But then the TF inferred distances and the mean 
peculiar velocities are subject to Malmquist bias. Suitable 
procedures to consider these biases, induced both by inho- 
mogeneities and selection function, have been discussed (see 
for instance Freudling et al., 1995, and references therein) 
where the spatial distribution, selection effects and observa- 
tional uncertainties are realistically modeled through Monte- 
Carlo simulations. We have adopted in our analysis inverse 
TF distances referred to the Cosmic Microwave Background 
frame (Wilhck et al., 1995; Wilhck et al., 1996; Willick et 
al., 1997). Inverse TF distances overcome distance depen- 
dent selection bias (see for instance Teerikorpi et al., 1998), 
nevertheless we have also tested the results using forward 
TF distances, fully corrected for Inhomogeneous Malmquist 
Bias by Willick et al. (1997). 



3 OUTLINE OF THE ANALYSIS 

The different methods applied to infer the distance to a 
galaxy are subject to uncertainties due to observational er- 
rors as well as scatter and systematics of the galaxy parame- 
ters. Distances are derived from linear relations between ab- 
solute magnitudes M and physical properties independent 
of distance as for instance the circular velocity Vc in the 
TuUy-Fisher relation, or the central velocity dispersion in 
the Dn — a relation. Both rms scatter and posible shifts in 
the zero-point of the distance relation should be taken into 
account in studies of the peculiar velocity field considering 
their strong influence on the results (Padilla, Merchan & 
Lambas, 1998). 

Since peculiar velocities of galaxies are inferred from 
redshifts and independently estimated distances, the effect 
of a zero-point shift would be observed as a systematic mo- 
tion of a shell of galaxies proportional to distance. The ef- 
fects induced in the velocity field by the scatter in the dis- 
tance relations (a or) depend on the catalog radial gradient 
which is affected by distance uncertainties. 

We correct the observed radial gradient by consider- 
ing a gaussian distribution of distance uncertainties. There- 
fore the resulting distribution of distance measurements of 
galaxies restricted to the same true distance bin da is ap- 
proximately gaussian centered at da-, with scatter ~ da (Jdr- 
Here da is given in units of km/s and (jdr corresponds to 
a distance fraction. Galaxies from other distances db [uab) 
will also contribute to the measured number of galaxies at 
da- The corresponding contribution from objects at dt can 
be expresed as: 



riab = T{db — da,a)nb, 



T{db-da,a) = 



(1) 



2Tva 



(d6-d„)+A/2 



e 2<T^ dx. 



where a = db aoR and A = 300km/s is the adopted bin- 
ning of galaxy distances corresponding to shells. Then, the 
number of galaxies measured at distance da (n^) takes into 
account contributions from all other distances: 



n'a = 2_^ ^"f' 

dh=0 



(2) 



where dmax is the limiting distance imposed to the catalog. 
In order to solve equation gI we define the vectors A^' = 
n'a and N — Ub and the matrix A = Uab- Then equation g 
can be rewritten as 



N' = AN 



(3) 



A can be inverted to obtain the true number count of galax- 
ies, A^, unaffected by the distance estimator scatter taken 
into account in the matrix A. However, a direct inversion 
of the matrix (using Gauss method for instance) produces 
diverging solutions for the last components of A*'. This diver- 
gence is produced by accumulation of large errors through 
the calculation over matrix rows. To avoid this problem, we 
used an iterative method in wich we apply A to A'^' and ob- 
tain N". Since A ~ /, the vector Ni = N' + {N' - N"), wiU 
be a better approximation to N than A^'. After k iterations, 
we obtain 



Nk+i =Nk + {N' - Ni) 



(4) 



where we impose the condition Nk > 0. 

This method also accumulates errors in the last com- 
ponents of A'fe, but the solutions start to diverge only when 
more than 3 iterations are applied. The optimal number of 
iterations was found to be fe = 2 or 3. The results of the iter- 
ative method have also been checked with the mock catalogs 
analysed in section 5, and were found to be accurate. 

In order to study the effects induced by the TF scat- 
ter in the velocity field, we calculate the effect of aoR 
on the mean radial velocity of galaxies {va) at distances 
d £ [da — A,da + A] under the assumption that the shells do 
not expand nor contract. We may write the apparent mean 
velocity of the shell v'a at distance da as: 



1 






(5) 



where v'^ is the individual velocity of the galaxy i. We recall 
the fact that at distance da there are contributions from 
other distances. If a galaxy j at real distance db is measured 
to be at distance da, the inferred velocity will be Vj — Vj — 
(db — da) where Vj is the real peculiar velocity of the galaxy 
j. Finally, if we sort by real distance galaxies accidentally in 
the shell at distance da, we can rewrite the last expression 
as: 



djnax ^ab 



'1 = -T ^ ^(uj - (db - da)) 



(6) 



db=0 J = l 

If we consider our assumption Vb — and add the possible 
presence of a zero-point shift Po, we find the final expression: 

(7) 



1 

v'a = —r / nab{da - db + Podb). 



(dt,-da)~A/2 



The inputs of this equation are the number count of 
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Table 1. Observations: The Mark III spirals 
Subsample 



of Gx. 




TF relation 


(JTF 


359 


Mh 


= -5.95 + 10.29r? 


0.47 


1355 


Mi 


= -5.79 + 6.8r? 


0.43 


383 


Mr 


= -4.28 + 7.12»7 


0.38 


156 


Mr 


= -4.18 + 7.73r; 


0.38 


326 


Mr 


= -4.22 + 7.73r; 


0.38 


433 


Ml 


= -5.48 + 7.87r; 


0.4 



Aaronson et al. Field (1982) 

Mathewson et al. (1992) 

Willick, Perseus Pisces (1991) 

Willick, Cluster Galaxy (1991) 

Courteau-Faber (1993) 

Han-Mould et al., CI. Gx. (1992) 




2000 4000 6000 8000 
d [km/s] 




2000 4000 6000 8000 
d [km/s] 



Figure 1. Figure la: Mean radial velocities of shells of galaxies from Mark III using inverse TF distances (solid lines) with a = 0.41 and 
Po = 0.05. The dashed line corresponds to the model results with the same values of a and Pq. Figure lb:Mean radial velocities using 
forward TF distances corrected for inhomogeneous Malmquist Bias (solid lines). The dashed line corresponds to the model results with 
parameter cttf = 0.3 and Pq = —0.13, and the dotted line to a model with cttf = 0.4 and Pq = 0. 



galaxies as a function of distance, the scatter aoR, and the 
zero-point shift, Pq- By comparing the measured values v'a 
from a catalog with the calculated values given by equation 
m we may infer the uncertainties that affected the measured 
distances. 

A similar deduction can be applied to obtain the appar- 
ent root mean square velocities a'^^ corresponding to galax- 
ies in a given shell at distance da- The assumption made 
here is that the true root mean square velocity of a shell is 
independent of distance. 

This quantity can be calculated from the following equation: 



^Idim 



Ed-mn^ / I I 2 X ^"^TTIUT r 7 7 . 7-7 7 l2 

rf^^o y^o-^^a - l^d^^O ^o.b[db -da + Podb\ 



dZZ < 



■(8) 



The inputs of eq. H are the observed o"^^, n^, the distance 
relation scatter odr, and its zero-point offset, Pq- The 1- 
dimensional velocity dispersion, (Tidim, is thus directly ob- 
tained from radial velocity data. The 3-dimensional velocity 
dispersion a^ is simply rj^ = *\/3o-idim, assuming isotropy. 



4 APPLICATION TO THE MARK III 
CATALOG 

The scatter of the Marklll spiral TF relation (hereafter (ttf) 
has been extensively studied (Mo et al., 1997; Willick, 1991; 
Mathewson, Ford and Buchorn, 1992), and similarly uncer- 
tainties of the TF zero-point, Po (Shanks, 1997; Willick, 
1991). The value of otf for the different spiral galaxy sam- 
ples in the Marklll catalog is around 0.4 measured in units 
of absolute magnitude (Willick, 1991). These authors give a 
null shift in the distance estimations with a mean deviation 
of ±0.07 in units of absolute magnitude. 

Since the sample of Mark III spirals has a reasonably 
smooth sky coverage the analysis outlined in the previous 
section is suitable for statistical purposes. We applied eq. 
R to the Marklll spirals restricted to distances d < 6000 
km/s since beyond this distance the uncertainties make pe- 
culiar velocities unreliable. We have calculated a x^ devia- 
tion between predicted and observed mean peculiar velocity 
of shells v'a as a function of axF and Pq. The best-fitting 
was obtained for 

axF = 0.41, Po = -0.05 

both in units of absolute magnitude. Notice that these values 
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Figure 2. Occurrence of values Pq and otf by random resamplings of the sample. Figure 2a: using inverse TF distances, Figure 2b: 
using forward corrected TF distances. 



of dTF and Po are consistent with those quoted in Willick 
et al. (1997). Figure hla and hlb show the observed mean ve- 
locities and the results of eq. ]A for the values of otf and 
Po obtained using inverse TF and forward corrected TF re- 
spectively. In figure lb it is also shown the predicted v'^ with 
no zero-point shift Po and otf = 0.4. In order to measure 
the accuracy of the determination of the distance uncertain- 
ties, we apply a x^ test to obtain otf and Po for a large 
set of catalogs obtained through bootstrap resampling. The 
number of realisations with resulting otf and Pq, are shown 
in figures Sa and Ep for inverse and forward corrected TF 
distances respectively. 

These results obtained for the inverse TF distances may 
be compared with the homogeneous calibration of the TuUy- 
Fisher relation corresponding to different samples of spirals 
of the Mark III catalog given by Wilhck et al., 1995; Wilhck 
et al., 1996; and Willick et al., 1997; which corresponds to 
values of otf in the range 0.38 — 0.47 mag (see Table 1). 

We have obtained the galaxy root mean square velocity 
by application of equation H to the sample d < QOOOkm/s. 
The value obtained is a„ — (660 ± 50) km/s, where the er- 
ror was obtained from 1000 bootstrap resamplings (Barrow, 
Bhavsar & Sonoda, 1984). The distribution of CTi, derived 
from the different bootstrap resamplings is shown in figure 

i 



5 TESTING THE METHOD WITH MOCK 
CATALOGS 

We test the results of our analysis using mock catalogs de- 
rived from the numerical simulation. We adopted a S7 = 0.5, 
Qa = COBE normalised CDM model which reasonably 
reproduces several statistical tests of large-scale structure 
such as cluster abundances, correlation functions, etc. This 
particular model requires no strong bias, so each particle of 
the simulation corresponds to a galaxy. 



We have considered 1000 random observers in the nu- 
merical simulation by defining cones with different positions 
and orientations in our computational volume. We have in- 
cluded the observed strong radial gradient wich corresponds 
approximately to a selection bias due to a magnitude limit 
cutoff in the data to the mock catalogs. This can be seen 
from the observed distribution of absolute magnitudes which 
is nearly gaussian with mean ~ M* (the knee of the Lumi- 
nosity Function) and a ~ 1.5mag. Nevertheless, for our sta- 
tistical purposes it is not necessary to adopt a Monte-Carlo 
model using the Galaxy Luminosity Function in the simu- 
lations. It suffices to reproduce the observed radial gradient 
in the numerical models through a Monte-Carlo rejection 
algorithm. Furthermore, we restrict the resulting number of 
particles of the mock catalogs to be equal to the number of 
galaxies in the observational sample. 

Observational errors in galaxy distance estimates were 
considered assuming gaussian errors in the galaxy absolute 
magnitudes of the TF relation with dispersion axF, corre- 
sponding to relative errors in distances -j- ~ ^^^. Namely, 
we assign to each particle in the mock catalog a new distance 
dneui = d{l + s) where s is taken from a gaussian distribu- 
tion with dispersion corresponding to the TF uncertainty. 
Then, as the particle peculiar velocities are inferred from 
the galaxy redshift and distance, v„sm = v — d s. 

The N-body numerical simulation was performed us- 
ing the Adaptative Particle-Particle Particle-Mesh (AP3M) 
code developed by Couchman (1991). The initial condition 
was generated using the Zeldovich approximation and corre- 
sponds to the adiabatic CDM power spectrum with Q = 0.5 
and JIa ~ 0. We have adopted the analytic fit to the CDM 
power spectrum given by Sugiyama (1995): 



P(fc) oc 



k /;n(H-2.34g) 
A \ 2.34(j 



(9) 



where A = [+3.89g -I- (16.1g)^ + (5.46g)^ + (6.71g)'']^/^ 
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Figure 3. Distribution of ay for the observational sample corre- 
sponding to bootstrap resampling. 



Figure 5. occurrence of values Pq and cttf in random resam- 
plings of a mock catalog with crj-p = 0.41 and Pq = —0.05 




r =770-1 17 



dotted lines 

CT =773±145 




(T„ [kni/s 



Figure 4. Mean radial velocities of shells according to models 
(solid line) and mock catalogs (dashed lines). Error bars corre- 
spond to fluctuations arising from different observers in the sim- 
ulations. 



Figure 6. Probability distribution of cr^ corresponding to mock 
catalogs with axp = 0.3 and Pq = (solid lines), crj-p = 0.41 
and Pq = —0.05 (dashed lines), and aj-p = 0.4 and Pq = —0.2 
(doted lines). 



k 

rh 



Mpc, r = 



is the mi- 



h exp{-ilB-^h/0.5nB/i1) ' 

crowave background radiation temperature in units of 2.7K, 
and Jls = 0.0125/i~ is the value of the baryon density pa- 
rameter given by nucleosynthesis. The normalization of the 
CDM power spectrum is imposed by COBE measurements 
using the value of as and h (the Hubble constant in units of 
100 km/ s /Mpc) chosen from Table 1 of Gorski et al. (1995) 
corresponding to an age of the universe to — 12 Gyr. The 
computational volume is a periodic cube of side length 300 
Mpc. We have followed the evolution of A'^ = 5 x 10* par- 



ticles with a 64^ grid and a maximum level of refinements 
of 4. The resulting mass per particle is 2.05 x 10^'^/i~^Mq. 
The initial condition corresponds to redshift 2 = 10 and the 
evolution was followed using 1000 steps. At the final step 
{z = 0) the linear extrapolated value of erg is compatible 
with the normalization imposed by observed temperature 
fluctuations in the cosmic background. 

As a test of the ability of equation m to reproduce the 
distorted values of mean velocity, we plot in figure W the 
resulting v'^ for the average of 1000 mock catalogs, each with 
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Figure 7. Probability contours corresponding to the observational sample d < 6000fcm/s. Dashed lines corresponding to diflFerent values 
of /3. The plus sign indicates the best values of a and Pq. 



distance uncertaintes included. The error bars correspond to 
the rms deviation of v' from the different observers in the 
numerical simulation. We adopt the values of arr and Po 
used in the construction of the mock catalogs to perform 
the calculation of v'^ through equation M. It can be seen 
in the figure the excellent agreement between the results of 
the calculation and the mock catalogs. We have applied a x 
method to derive otf and Po- We find that our method gives 
accurate estimates of the actual arr and Po when applied 
to the mock catalogs. In figure H we show the frequency 
of aTF and Po obtained from bootstrap resampling objects 
corresponding to a mock catalog with an imposed scatter 
(jtf = 0.4 and Po = —0.05 as derived from the observations. 

The values of otf and Po inferred from the x^ method 
can be used to calculate the value of (t„ in the mock cat- 
alog through equation w. The distributions of inferred (t„ 
found are shown in figure H where it can be appreciated the 
stability of the method in estimating aidim under different 
zero-point offsets Po and (jtf- The actual value of the rms 
peculiar velocity of the simulation is (t„ ~ 780 km/ s which 
may be compared to the inferred values from equation H for 



different adopted values of cttf and Pq. For instance, we 
find a^ = (743 ± 102) km/s, (t„ = (770 ± 117) km/s, and 
cr„ = (772 ± 145) km/s for gtf = 0.3, Po = 0; gtf = 0.41, 
Po = —0.05; and axF ~ 0.4, Pq — —0.2 respectively indicat- 
ing the ability of our procedure to derive the rms peculiar 
velocity irrespectively of TF scatter and zero-point offset. 

We find a zero-point shift Po — —0.05 in our analysis of 
the observations using inverse TF distances. It is of interest 
to test the probability of occurrence of such a value arising 
from our model assumptions such as Va ~ 0, etc. We use 
mock catalogs with imposed a = 0.4 and Po = in the TF 
calibration and we test the probability of finding different 
values of zero-point shifts . We apply a x^ method to derive 
the pair of values gtf and Po that provides the best-fitting 
of equation m to the actual values for each mock catalog. We 
compute the frequency of occurrence of (jtf and Po for the 
difference observers, and we can estimate the probability of 
obtaining different values of zero-point shifts. We find that 
a random occurrence of the observed value Po = —0.05 is 
within a standard deviation, consistent with Willick et al. 
(1997) estimate. 
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6 DETERMINATION OF THE p PARAMETER 

The 3-dimensional velocity dispersion a^ is directly related 



to th e quantity P — Q ' /b through the relation (Peebles 



198C) 



-{¥)' 



y^{y)dy 



(10) 



where H is the Hubble constant, a is the expansion factor of 
the universe, y is expressed in Mpc, ^{y) the galaxy spatial 
correlation function, / ~ n°'® the rate of growth of inho- 
mogeneities, and b is the linear bias factor. We express y in 
units of Mpc h~^ , therefore H = WOkm/s/Mpc and a = 1. 
We estimate Q''^'^'/b from the inferred value of 3- 
dimensional root mean square peculiar velocity. We adopt 
the power-law fit to the galaxy spatial correlation function 
estimated by Ratcliffe et al. (1997): 



^(r) = 



(-)" 



if r < 50 Mpc 
if r > 50 Mpc 



with ro — 5.1 Mpc and 7 = 1.6. 

We calculate (3 from equation hd using equation M to 
express a^ and therefore l3 in terms of the parameters arp 
and Po- In figure mwe show equal /3 contours in the a ~ Po 
plane for our sample of Mark III spirals with distances d < 
6000 km/ s. Also shown in this figure are the 1 a and 2 
(T contour levels corresponding to the frequency of inferred 
c^TF and Po from bootstrap resamplings of the observational 
data set. The corresponding result is fi — 0.601qq5. 

The true uncertainty of the global value /3 — ff'^/b 
derived from a catalog of peculiar velocities will be grater 
than that obtained from bootstrap resampling of the data 
due to cosmic variance. A suitable value of the uncertainty 
in /3 can be estimated from the rms values of this parameter 
derived from the mock catalogs. According to our analysis 
A/3 ~ 0.1 for a limiting distance dmax < 6000. Thus, adding 
in quadrature both errors, we find /? — 0.601qj^^ for d < 
6000fcm/s. 

Thus, the observed peculiar velocity field is inconsistent 
with a critical density universe if optical galaxies trace the 
mass. 



7 CONCLUSIONS 

We have developed a method for the analysis of the peculiar 
velocity field inferred from peculiar velocity data and we 
apply this procedure to the spirals of the Mark III catalog. 
We estimate optimal values of inverse TuUy-Fisher scatter 
and zero-point offset for a sample of the catalog with limiting 
distance diim = 6000fem/s. We derive the 3-dimensional 
rms peculiar velocity of the galaxies a^ — (660 ± 50)km/s 
where the uncertainty has been obtained through bootstrap 
resampling of the data. 

In our model, the shells have not a net mean radial mo- 
tion, and the mean square velocities of galaxies in difi'erent 
shells are described by a unique number cr„. The compari- 
son with mock catalogs derived from numerical simulations 
shows that these are reasonable hypotheses that allow to ob- 
tain physical characterizations of the nearby universe from 
observations. We have shown that corrected TF distances 
require a large Po — — 0.15mag so that caution should be 
taken when they are used in statistical analysis. 



We use mock catalogs derived from numerical simula- 
tions of CDM models considering measurement uncertain- 
ties and sampling variations to check our statistical analysis. 
We find a general good agreement between the results of the 
calculations and those measured in the mock catalogs. The 
spread of a„ measurements from different observers in the 
numerical simulations may be added in quadrature to the 
bootstrap resampling errors to provide a more reliable esti- 
mate of the uncertainty in (t„. We infer a^ = 660 ±70 km/ a, 
and we conclude that (3 ~ OP'^ /b = 0.60to]_l- 

Estimates of the parameter /3 from other analysis such 
as studies of redshift space distortions of the galaxy two 
point correlation function provide similar values /3 ~ 0.5 
(Ratcliffe et al. 1997). Moreover, the confrontation of ob- 
served cluster abundances with prediction of different cos- 
mological models put constraints of the form Q" /b ~ 0.4 — 
0.6, values consistent with our determinations. 
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